Dipole Excitation of Dipositronium 
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The energy interval between the ground and the P-wave excited states of the recently discovered 
positronium molecule Ps2 is evaluated, including the relativistic and the leading logarithmic radiative 
corrections, E P — E s = 0.181 586 7(8) a.u.. The P-state, decaying usually via annihilation, is found 
to decay into the ground state by an electric dipole transition 19 percent of the time. Anticipated 
observation of this transition will provide insight into this exotic system. 
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Last year's discovery [l,] of dipositronium (PS2) was 
welcomed as a herald of a new kind of chemistry. Ps2 is a 
bound state consisting of two electrons and two positrons. 
Its stability against dissociation was predicted in Q| , but 
its observation was very challenging. Ps2 rapidly annihi- 
lates producing photons similar to those from atomic Ps 
decays. Nevertheless, the evidence of the Ps2 existence 
is now compelling. 

In the experiment described in [H, an intense pulse of 
positrons is stopped in porous silica, forming Ps atoms, 
some of which make their way into the voids of the pores. 
Ps atoms have two hyperfine states: a short-lived spin- 
singlet para-positronium (pPs) with a lifetime in vacuum 
of Tpp s = 0.125 ns, and a long-lived spin-triplet ortho- 
positronium (oPs), with r p s = 142 ns. Interactions may 
shorten the oPs lifetime through two mechanisms: spin 
exchange quenching (SEQ), in which spins flip convert- 
ing oPs into the rapidly-decaying pPs; and formation of 
molecules Ps2. In the latter case, as we will discuss below, 
the probability of each of the electron-positron pairs to be 
a spin-singlet is one quarter, and the size of the molecule 
is similar to that of an atom. Thus Ps2 is short-lived, 
with the lifetime of about 2r p p s . 

But even if the rapid disappearance of oPs is observed, 
how can one tell whether this is due to molecule forma- 
tion rather than SEQ? The key is that the molecule for- 
mation needs another body, such as the pore surface, to 
absorb the released binding energy. As the temperature 
of silica is increased, the fraction of Ps atoms on the sur- 
face decreases, fewer molecules should be formed, and 
more oPs survive. Exactly such an effect is observed []]]. 
More recently, evidence of the Ps2 formation on a metal 
surface has also been found 0. 

As well as proving the existence of the first known sys- 
tem containing more than one positron, this discovery is 
viewed as an important step towards studies of even more 
exotic phenomena: larger polyelectrons, Bose condensa- 
tion of positronia, and eventually a 7-ray laser based on 
stimulated annihilation. 

Before those exciting possibilities are explored, a more 
detailed study of the newly-discovered PS2 is warranted. 
Like the relatively better-known positronium ion Ps - 
[1, 03], the molecule PS2 is weakly bound 0, 0|- How- 
ever, whereas the ion has no excited states stable against 
dissociation, the molecule has an interesting spectrum of 



three: two excited S states and one P 0, 0]. The latter 
has a sizable branching ratio for an electric dipole transi- 
tion to the ground state, accessible to middle-ultraviolet 
laser spectroscopy [l(J. This is a great asset since poly- 
electrons generally decay through annihilation into pho- 
tons, too much like atomic pPs to reveal the structure of 
the decaying system. Anticipated observation and precise 
measurement of this line [llj will confirm the presence of 
the molecules and test our understanding of their nature. 

Measurements of the P-state excitation energy are ex- 
pected to have a precision of 10 parts per million (ppm) 
sensitive to relativistic corrections. The relativis- 
tic effects have been found recently for the ground state 
[12] but the P-level energy is only known in the non- 
relativistic approximation [8J]. 

The challenge of precise theoretical studies of PS2 is 
that it is a relatively complicated four-body system whose 
components have equal masses. Unlike in molecules built 
of atoms with hadronic nuclei, no part of PS2 can be 
treated as slowly-varying and the Born-Oppenheimer ap- 
proximation cannot be applied. Nevertheless, its ground 
state has already been thoroughly studied with varia- 
tional and quantum Monte Carlo methods (see for a 
recent review of earlier work) . The P-state is more chal- 
lenging. Because its wave function has a node, the varia- 
tional procedure converges slower. Larger expressions for 
the matrix elements and additional integrals exacerbate 
the difficulties [3, [Hi . 

To overcome these obstacles, we combined the varia- 
tional method in a Gaussian basis with algorithms for 
decomposing the Hamiltonian matrix and for optimizing 
the wave function. In addition, we sped up the conver- 
gence by transforming the operators representing correc- 
tions to the energy, using a method proposed by Drach- 
man As a result we not only match or exceed the 

accuracy of the best existing evaluation of the relativistic 
corrections to the ground state, but also extend those re- 
sults to obtain the leading next-order corrections (QED) 
and, more important, determine analogous effects for the 
P-state. We find the energy interval 



AE = E P - E s = 0.181 586 7 (8) a.u., 



(1) 



or 4.941 23 (2) eV, corresponding to the wavelength A = 
250.9179(11) nm. The branching ratio for the dipole 



2 



transition is also determined, 

r dip (P -> S) 



BR(P-> S) = 



Tannih (P) + r dip (P - S) 

The dipositronium's Hamiltonian is 

H = He + a 2 H Te} + O (a 3 In a) 



0.191(2). 



(2) 
(3) 



Its leading term describes the non-relativistic Coulomb 

system, 



He 



a<b 



Zab 
Tab 



(4) 



where p a and r a b = \r a — f%\ are momenta and relative 
distances of positrons (a, b = 1, 2), and electrons (3, 4). 
z a b equals +1 for a like-charged pair a, b and —1 for oppo- 
site charges. As units of length, momentum, and energy, 
we use l/am e , am e , and a 2 m e , and set c = K = 1. 

Wave functions and the lowest-order ("nonrelativis- 
tic") values of energy levels are determined by the 
Coulomb Hamiltonian |(4|). Higher-order corrections 
O (a 2 , a 3 In a) to energies are computed as first-order 
perturbations with those wave functions. 

The Hamiltonian J3]) has a rich symmetry 0, Q3] that 
is reflected in the wave functions. In addition to the sym- 
metry with respect to permuting coordinates of electrons, 
as well as those of positrons, there is also the charge con- 
jugation symmetry. If the system is described by relative 
coordinates r a b only, the latter is equivalent to the spa- 
tial inversion, and described in a given state by its parity 
7T. Thus the wave functions of the ground state with 
L 7 * = + and of the P-state 1~ are ipi — A\x4>i], where 
X = | (TiJ.2 - I1T2) (hU - I3T4) is constructed using 
electron and positron spinors and the antisymmetrizer 
is built out of operators permuting pairs of particles, 
A = 4? (1 + ttP 13 P 2 ±) (1 - P l2 ) (1 - P 34 ). The spatial 
are expressed in a Gaussian basis, 



1 

wave functions 



iV 



1=1 

JV 



a<b 



t>p = ri Ci exp 

i=l 



b 



a<b 



Here fi denotes a particle coordinate with respect to the 
center of mass. Since r*i + f 2 + ^3 + ^4 = 0, in terms of 
the relative coordinates we have r\ = -r (f\% + ri3 + F14), 
The six parameters w^' P in are determined, for 
each of the N elements of the basis, in an extensive op- 
timization process. QR decomposition [l3| and inverse 
iteration are used to determine energy eigenvalues of He- 
In the i th step of minimizing the energy, the six parame- 
ters of the i th basis element are optimized with Powell's 
method without gradient [l8j]. The optimization steps 



are cycled through the basis elements until convergence 
is reached. This procedure yields nonrelativistic energies 
accurate to better than one part per billion, shown in the 
first line of Table U The ground state agrees with 
although our error bar is slightly larger. For the P-state, 
we improve on the previous best result @|. 

To test the numerical procedure, we used the ground- 
state of lithium, presently the best known four-body sys- 
tem. Its ground-state energy was evaluated using Hyller- 
aas coordinates fTsl . lit! l2l| with a precision of about 
10~ 12 . It reliably calibrates both the absolute value and 
its uncertainty since, thanks to the richer symmetry of 
PS2, the Gaussian method converges better for it than 
for Li. 

Assured of the quality of the obtained wave function, 
we proceed to the relativistic corrections. For the two 
states of interest, the following parts of the Breit-Pauli 
Hamiltonian [3 contribute, 



ffrel = H M V + H D + Hqo + H SS + H A , 

a 

Hu = -Try^ z ab S 3 (rab) , 



a<b 



Hoo — 
Hss = 
H A = 



1 a<b V ab ) 



pI- 



2- 



^ Zab d a ■ S h S 3 (rab) , 



(6) 
(7) 

(8) 

(9) 
(10) 



a<b 



| ]T (3 + a a -a b )d 3 (r ab ). (11) 



a<6,afc^l2,34 



Given the spin configuration of PS2, S a -a b can be replaced 
by zero for an e + e~ pair, and by —3 for like-charged 
pairs fl3 ] . Thus, only four operators remain: p 4 , delta- 
functions for e + e~ and e~ e~ , and -ffoo- 

A disadvantage of the Gauss basis is its incorrect be- 
havior at short inter-particle distances: it does not re- 
produce the cusps of the wave function. This slows 
down the convergence of matrix elements of the delta- 
function and the kinetic energy ("mass-velocity") cor- 
rection p 4 . For example, in lithium with a basis size 
(5) of 2000, the error is a few units in 10 4 Q. Even 
more dangerously, the convergence is so slow that a mis- 
leading limit may be deduced. To overcome this diffi- 
culty, the operators can be transformed into an equiva- 
lent form, whose behavior is less sensitive to the short- 
est distances. For the delta-function, a prescription was 
found by Drachman Neglecting boundary terms, 

^S 3 (r ab ) <h<h - £ (E - V) frfo + E c (V ^r^ 2) , 
where V = ^2 a <b~- For the kinetic energy we use 
J2 a fa^fa - 4 (E -V) 2 E a<fc (V 2 Ji) (V 2 2 ) . 

Numerical results for the four basic relativistic opera- 
tors are shown in Table HH For the ground state we find 
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agreement with [12| to within 1 ppm for i?oo, and three 
or four digits for the remaining operators. Since these are 
the operators that we regularized, we believe our results 
to be more accurate, despite the much smaller size of the 
basis used (the difference is unimportant since unknown 
higher-order effects are likely larger). Again, lithium was 
used as a test. The kinetic correction converges slowest, 
as usual in many-body calculations [13]. 

Table U shows corrections to the energy levels. The 
orbit-orbit and the Darwin terms largely cancel. The 
result is dominated by the virtual annihilation and, to 
a smaller degree, by the kinetic term. The annihilation 
is repulsive and decreases the binding. It is about twice 
as effective in the ground state as in the P-state so that 
overall the relativistic effects decrease the S-P interval. 

Beyond relativistic corrections, the next largest effect 
is of the relative size O (a 3 In a) , analogous to the Lamb 
shift in hydrogen. Its physics is richer in positronium 
[i3 ] because two-photon "recoil" interactions between two 
light constituents contribute at the same order as the 
self-interaction corrections. In hydrogen the latter dom- 
inate, the recoil effects being suppressed by the electron- 
to-proton mass ratio. In Ps2, one should in principle con- 
sider interactions among all pairs of particles. However, 
contributions of like-charged particles are suppressed by 
two orders of magnitude, as can be seen by comparing 
the last two columns of Table Ull Thus the coefficient of 
a 3 In a is given with excellent accuracy by expectation 
values of a 3 In a i/i og = —24a 3 In a 5 3 (ris) [25(| , shown 
in Table HI We take halves of their values to estimate the 
error for each level, and add those to obtain the error 
estimate of the final result, Eq. (!]). 

The precise value of the e + e~ overlap given in Table ITT1 
provides a new prediction of the annihilation rate of PS2 
in both states, 



the annihilation. The rate of the dipole transition is 



r a „nih(S) = 4™ 3 (<5 3 (r 13 ))(l + RC) 

= 1/(0.22455(6) ns), 
r annih (P) = 1/(0.442 77(11) ns), 

where the RC denotes the radiative corrections, 



(12) 
(13) 



RC = a 



19?r 
~12 



-)+2a 2 ln-+0(a 2 ) 
it J a 



The error, due to unknown O (a 2 ) effects, is estimated 
as half of the logarithmic correction. The annihilation 
in state S was previously calculated in [2^]. Our result 
differs slightly, primarily because of the error in the O (a) 
correction in that study. The result 11121) confirms the 
expectation that the lifetime of PS2 in its ground state 
is about half of that in pPs. In the P-state, one of the 
particles overlaps only negligibly with the others, slowing 
down the annihilation by another factor of two. 

For the experimental search of the dipole transition, it 
is interesting to know how competitive it is relative to 



r dip (P -> s) 



-a 3 (E F 



1/(1.873 ns) , 



(14) 



where d = f*i + — rs — r±. The dipole matrix ele- 
ment is determined as a Gaussian integral, (s d p\ = 

2.040 942 265(16), and leads to the final result for the 
branching ratio, Eq. Corrections to this prediction, 
conservatively estimated as less than one percent, arise 
from the three-body decay into a photon and two Ps 
atoms, and from relativistic effects O (a 2 ) [13]. The 
present value exceeds the previous evaluation [l3] by 
about ten percent, a welcome improvement for the exper- 
imental search of this transition. The larger value may 
slightly facilitate the use of the P-S transition intensity 
to determine the number of the PS2 molecules produced. 

It is noteworthy that the rate of the dipole transition 
is similar to twice that in a positronium atom. This con- 
firms the approximate picture of the excited P-state as re- 
sembling two weakly interacting positronium atoms, one 
in its ground state and the other in the P-state Q. In a 
molecule consisting of two weakly interacting atoms, the 
P-S dipole matrix element is, because of coherence, y/2 
times larger than in an isolated atom. This approxima- 
tion predicts |(cZ)| = = 2.1. Similarly, in molecular 
hydrogen one finds [2I \{d)\ = fj§ = 1.05, in the limit of 
weak interaction between the two atoms. 

Fortunately, the 2P-1S energy interval does differ suf- 
ficiently between atomic and molecular positronia for 
its measurement to unambiguously confirm the existence 
of PS2. The main difference arises already in the non- 
relativistic energy values. In the molecule, one Ps atom 
may be interpreted as a dielectric medium that weak- 
ens the electric field in the other one, thus decreasing all 
energy intervals. Relativistic effects, primarily the anni- 
hilation, slightly add to that decrease. The dipole matrix 
element is also decreased below the asymptotic value of 
V2 times that in a free atom, weakening the transition 
rate below half of the atomic rate. 

On the technical side, this study reveals the somewhat 
unexpected potential of the correlated Gaussian basis. 
The fast optimization method described here leads to 
comparable or better results than previously published, 
even with a much smaller basis. A drawback of Gaussians 
is their incorrect asymptotic behavior, both at short and 
at long distances. This is compensated by the availability 
of an analytical form of all required matrix elements and 
by the good numerical behavior of the integrals. Double 
precision sufficed for the variational parameters. 

The Gaussian basis is especially suitable for the 
positronium molecule since it tracks all inter-particle dis- 
tances. The high symmetry of PS2 improves the conver- 
gence of the variational procedure. Parallelizing the code 
would certainly lead to improvements, and will likely 
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Source 


Ground state 


P state 


H c 


-0.516 003 790 415 (86 


) -0.334 408 317 34(81) 


a 2 HMV 


-0.000 009152 


-0.000 004 780(1) 


a 2 Hoo 


-0.000 013 470 


-0.000 007 736 


a 2 H 


0.000 014 592 


0.000 007458 


a 2 Hss 


0.000 000 419 


0.000 000 097 


a 2 Ha 


0.000 022 202 


0.000 011259 


2 TT 

a H ve \ 


0.000 014 591 


0.000 006 298(1) 


a 3 In a H\ og 


0.000 00101(50) 


0.000 000 51(25) 


Total 


-0.515 988 2(5) 


-0.334 4015(3) 


[12] 


-0.515 989199 656 





TABLE I: Corrections to the energy levels of Ps2. 



be necessary for the determination of further QED cor- 
rections. For the present and foreseeable measurement 
goals, the theoretical description of the dipole transition 
energy and its probability presented here is sufficient. Its 
experimental test will complement the newest chapter in 
chemistry with one in spectroscopy. 
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Pachucki, Alexander Penin, and Mikhail Voloshin for 
very helpful discussions. We thank Alexander Brown 
and Lutoslaw Wolniewicz for advice on the molecular hy- 
drogen literature, and Paul McGrath for improving the 
manuscript. This research was supported by Science and 
Engineering Research Canada. 
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Basis size 


<E„i*> 




10 2 (S 3 (r u )) 


10 4 (<5 3 (r 12 )) 


Ground state 


2200 


1.374 923(45) 


0.505 892 400(27) 


2.21185117(14) 


6.256 827 3(42) 


1600 [10] 






2.21151 


6.259 


6000 [121 


1.374 696 3 


0.505 892 40 


2.211 775 9 


6.257 950 5 


State P 


2200 


0.718150(86) 


0.290 557 920(46) 


1.121 723 38(31) 


1.453 512 7(82) 


1600 






1.12091 


1.4591 



TABLE II: Expectation values of the basic relativistic operators, compared with previous studies, where available. 



